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1 Introduction 

The Skyrme model [2] was proposed to repair weak points of QCD at low energy (in GeV) . Time 
by time, it really become a phenomenological model of baryons. In 1990, H. Y. Cheung, F. Gursey 
[1] gave the general Skyrme model in which the sigma model [3] and the V — T model [9] were two 
special cases corresponding to n—1 and n—2. When Cheung and Gursey constructed composite 
Skyrme model, they also didn't let the mass of pion into the Lagrangian density [1] because of its 
spontaneous symmetry breaking [3, 4]. However, datum of prediction are not really axact. However, 
datum of prediction are not really axact [1]. So, we have proposed the Lagrangian density with the 
pion's mass. In 1984, E. A. Bergshoeff, R. I. Nepomachie, H. J. Schnitzer proposed Skyrmion of 
four-dimensional supersymmetric non-linear sigma model [5]. However, there were not papers that 
concerned with an extending the composite Skyrme model to the supersymmctry. Thus, we have 
performed to extend the commposite Skyrme model into supersymmetry. Recently, the Einstein- 
Skyrme system has been studied by several authors [6, 7, 8]. And then the skyrmion to be the 
gravitating skyrmion. All of these study were related to the non-linear Skyrme model but were 
not related to the composite Skyrme model. Thus, we have constructed the gravitating composite 
skyrmion. 

2 The mass of pion 

The Lagrangian density is defined as 

C*T = 4. + P> (U n ) - 2] . (1) 

Putting the hedgehog solution given by Skyrme [2] U a (r) = exp [ir.fF (r)] into the La- 
grangian density (1), we have 

= Cn + (cos nF - 1) , (2) 
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where r's are Pauli's matrices, f = ™ and 
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The static energy of chiral field is 
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and f = eF„-r is a dimensionlcss variable. 

Now, we shall quantize the hedgehog solution ( or soliton) by collective coordinates A (t) = 
a (i) + ia (t) f with coditions a^a^ = 1 and a^a^ = (u = 0, 1, 2, 3). By this way, solitons being 
real partilces. Different quantizations make different real particles. With the A (t) matrix, the 
U n (x, t) will transform as U n (x, t) = A (t) U$ (x) A- 1 (t). We put U n (x, t) into Eq. (1) then the 
Lagrangian is 

3 

L n = -£r n + T„Tr [d A (t) dvA- 1 (t)] = + 2T n ]T d 2 , 



where 



e 3 F w 
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sin 2 nF < 1 + 4 
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The conjugate momenta are defined as 



ju=0 



/ddf, = 4r„d AI . 



(10) 



From Eq. (7), we can get the Hamiltonian 
H n = tt^ - L n = AT n a^ -L n = £%° n + 2T n d^ = + ^ ^ /8T n . (11) 

Canonical quantizing momenta ir^ = —id/da^ takes the Hamiltonian as 



H n = £r n + (l/Wn)J2(-d 2 /dal). 



(12) 
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The isotopic vector c is defined as 

-A* A = (do — idjTj) (a + ia,Tj) = a a + aa — i ^a a — d a + a x ~aj f = — zcr. (13) 

The square isotopic vector c is 
c 2 = a 2 ciQ + a 2 eiQ + a^a 2 — (da) 2 — 2doaoaa = a 2 do + a 2 ao + aPa 2 — (da) 2 + 2(aa) 2 = d M d M . (14) 

In the configuration of hedgehog, the spin and isospin rotation are equivalent. We can 
homogenize c 2 with the square of spin and isospin 

* = ^ = I^ 2 - M 
Thus, we have the Hamiltonian 

Hn = + 7^ J 2 = £ P n° n + 7^T 2 , (16) 

81 „ 81 „ 

where J is a spin and T is an isospin. Eigenvalues of Hamiltonian (16) are 

M n = E n = ^A n + -^-Q n + ^±^, (17) 

where I = 2 J = 21 and (J, /) are respectively the isospin and spin quantum numbers. So, the mass 
of nucleon (J = 1/2) and delta (J = 3/2) are given by 

M " = f A " + ?k e " + W' (18) 
"-T^ + ^ + W- (19) 

From Eqs. (18) and (19), we can get 

2 (Ma - Miv) , . 
^3 > ( 20 ) 



4(M A -M N y$lA n 



6M N (M A - M N f $„ - 96„ - 3 (M A - M N f * r 



(21) 



We obtain the nonlinear differential equation of F (f) from the minimum condition of the 
chiral field's energy: 5 F £ P ion = 

/ f 2 2 \ ., f . sin 2nF / , sin 2 nf 1 \ „ 9 . , 

— + -j sin 2 nF F" + -F + F' 2 - - - /3r 2 sin nF = 0, 22 

\ 4 / 2 n \ n z r z 4 / 

where 
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We see that Eq. (23) contain the constant of e and but their values are not known. 
They can be found after we solve Eq. (22). Thus, we must use the special way to solve Eq. (22). 
Firstly, we assume the value of (3* , put it into Eq. (22) to solve this equation. When numerical 
solution of Eq. (22) is found, values of F v , e and (3** will be found from Eqs. (20), (21) and (23). 
After, we shall compare two (3* and (3** until they equal approximately together. For the case of 
n = 4, we need the value of (3 is 0.0025. 

After constructing coposite Skyrme model with pion's mass, we shall apply this model to 
calculate some static properties of nucleon (proton and neutron) [1] in a case of n = 4. And 
numerical results of static properties of nucleon are described in Fig. (2) and the below table 

NUMERICAL RESULTS TABLE IN A CASE OF N=4 

Quantity Prediction Experiment 

e 2.87 

186 
0.71 

0.88 

0.81 

0.8 
2.79 
-1.91 
1.46 
1.23 

3 The supersymmetric composite Skyrme model 

The Lagrangian density of composite Skyrme model is defined as 

Cn = St Tr (W ndtlUn ) + 32kl Tr (\- U ~ nd » Un > U ^ Un ] 2 ) > ( 24 ) 

where is the pi-meson decay constant, e is a dimensionlcss parameter. 

The ordinary derivatives in the Lagrangian density (24) is replaced into the covariant deriva- 
tives 

D„U n = d„U n - iV^U n r 3 , (25) 

then Eq. (24) is 

f 2 



i^(MeV) 


446.56 




0.21 
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^„ eu (mag) 
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1.22 
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Tr(D^U- n D^U n ) + ^L^Tr ([U- n D^U n ,U- n D u lI n ] 2 ) . (26) 
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or 



Eq. (26) is invariant under the local U (1) R and the global SU (2) L transformations 

U n (r) -» AJ7 n (x) e iA « r3 , A e St/ (2) L , (27) 
V»^(r)+^A(r). (28) 
The gauge field (r) in Eq. (28) is defined as 

V^ = -^Tr{U- n d„U n T 3 ). (29) 

We parametrize the SU (2) matrix U n (r) in terms of the complex scalars Aj 

^w=(l; - A f ). (so) 

From the unitary constraint, we have 

U^U = 1 => A l A l = A\Ai + A* 2 A 2 = 1. (31) 

Eq. (25) can be rewritten as 

A 1 -A%\ _ ( Ai -At, \ ( A\ A* 2 



D ^ 2 a{ =M£ At - iV »{2 -A ) (32) 



D li A i = (d lt -iV lt )A i , (33) 
D^A i = (d^ + iV^)A i . (34) 



Then, the form of gauge field (29) is 



i 



V^{r) = -—A*dA i . (35) 
In 

With the matrix U n (r) defined in Eq. (30) and the gauge field defined in Eq. (35), the 
supersymmetric Lagrangian density is 

where 

F\xv (V) = dfjYv — dvVfj,. (37) 

To supersymmetrise composite Skyrme model, we extend Ai to chiral scalar multiplcts 
(Ai,tpai, Fi) (i,a= 1,2) and the vector V^(r) to real vector multiplcts (V^,A Q ,D). Here, the 
fields Fi are complex scalars, D is real scalar, ip ai , A Q are Majorana two-component spinors. ip a i 
corresponds to a left-handed chiral spinor, tp m — (V>f )* corresponds to a right-handed one. The 
supersymmetric Lagrangian density is given by 
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r - Jx. 

J~susy — g„2 



-D^D^Ai - \ii>™ [p v ) a . D xjjf + F*F t - 



-iA l X a ^ m + iA^i + D {A l A t - 1)] + 



+ 



8e 2 n 2 



It is invariant under the following set of supersymmetric transformations [5] 

5A t = -e a tp m , 

5ip ai = -ie* (o- m ) qA D^A, + e a F t , 
5F> = -ie« [a^l - is^A'X^ 

5D= l -(a») a .d,{^X a ~e a X«). 



(38) 

(39) 

(40) 
(41) 

(42) 

(43) 

(44) 



The field equations and their supersymmetric transformations lead to the following con- 
straints [5] 

A'Ai = 0, 



A'lfiai = 0, 

A l F t = 

and following algebraic expressions [5] 

K„ = -^{a i 5 /1 ^ + (^) ad ^Vai}, 

\ a = -i {Fty aj + i K) QA (D^Ai) , 
D = D^AW^ + (p^fj - F l Fi. 

Setting ip a = F t = 0, Eq. (39) is 

1 



fl - - 1 

Zsusy ~ ~^D»AD^A + ^2 



(45) 
(46) 
(47) 

(48) 

(49) 
(50) 

(51) 



We see that in Eq. (51) the second term is fourth-order in time derivatives. However, there 
are other possible fourth-order terms [5] 



a ADA - (d^ad^a) 2 . 



Thus, we can add Eq. (52) into Lagrangian density (51) 



(52) 



■'Susy 



-iKD^AD u A 
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+/3{nAnA- (DM£>„i4) 2 }] . (53) 
The hedgehog ansatz is defined as [1, 2] 

T 

U n (r) = cos nf (r) + if- sin nf (r) . (54) 

From Eq. (30), we have 

A\ = cos nf (r) + i cos 6 sin nf (r) , (55) 
A 2 = ie llp sin 9 sin nf (r) . (56) 

Inserting Eqs. (55) and (56) into the supersymmetric Lagrangian density (53), we obtain 
the static energy 



+ ^(/" + ^-^]!. (57) 



From a minimum condition of the hedgehog's energy SfE — 0, we obtain the following 
nonlinear differential equation of f(x) 

-x 2 f" - 2xf + 2iB|!Lf + Hs+Sl \ 2 f"if n f - 6a; 2 /' 2 /" - 4x/' 3 + 



, /' sin 2nf , sin nf sin 2n/ 

n n 3 a: 2 



+ 



4/' 2 sin 2n/ 4 sin 2 n/ sin 2n/ 



x 2 f(4) , 4x/< 3 ' _ 4/"cos2n/ , 
•* 1 77. tj. 1 



0, (58) 



where x = e/^r is the dimensionless variable. 

Two boundary conditions of Eq. (58) are defined as 

/ (0) = tt, (59) 

/ (oo) = 0. (60) 

Eq. (58) can not give the analytic solutions. So we use numerical method to solve it. 

4 The gravitational composite Skyrme model 

The Lagrangian composite Skyrme model coupled with gravity can be defined as 

£ = £g+£cs (61) 
where Lg = 16 * G _R, R is a curvature scalar and G is a Newton constant, 

Us = -J^d^Tr (8^U n d v U n ) + (62) 
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+ m^^g^Tr [(d^U n ) U~ n , (d p U n ) U~ n ] [{d v U n ) U~ n , (d a U n ) U~ n ] , (63) 

U(x) is an SU(2) chiral field, is the pion decay constant, e is a dimesionless parameter. We 
consider the static spherically symmetric metric given by 

1 



ds 2 = -N 2 (r) C (r) dt 2 



C(r) 



dr 2 + r 2 (d0 2 + sin 2 9d V 2 



where 



From the equation 



we can define tensor metric g pv 



C(r) = l- 2 -^. 

r 

ds 2 = g^dx^dx 1 ' , 



(64) 

(65) 
(66) 



<fc„ = (-N 2 (r)C(r). 



,1,1 



diag 



C(r) 

Inserting the hedgehog ansatz proposed by Skyrme [1] U (r) — cost (r) + mrsinf (r) into the La- 
grangian density (1), we can obtain the static energy density for the chiral field 



(67) 
(68) 



£ s = -JL ( Cf' 2 + 2 



sin 2 n/\ 1 sin 2 nf f ori ,/2 sin 2 nf 



+ 



2c.r + 



(69) 



where r are Pauli matrices and n = - . 

Introducing dimensionless variables, they are useful in solving a field equation 

X — & 7j- T* . 

u (x) = eF n m (r) . 
In terms of x and fi (x) , the static energy can be written as 

Es = J y/^gSs^n^dr, 



or 



where 



r = 



e/' : 



E s = ^—T, 

e 



sin n/\ sin nf 



2n 2 x 2 



2Cf' 2 



sin nf 



Nx 2 dx. 



The energy tensor is defined as 



T- 



'-g 



where 



dx k dg l k 3 



(70) 
(71) 

(72) 
(73) 

(74) 
(75) 
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Inserting the form of Lagrangian density (61) into Eq. (76), we can get 



Ty = " 2^ F ' Tr ( LiLj ~ bn^lA + ^Tr [Fi^g^ - \ 9lJ F^ ) ( 77) 



where 

L„ = U- n d^U n , (78) 

F^v = [Lp, L v ] . (79) 
The Einstein's gravitational field equation is 

Ri j -^g ij R = 8nGT ij , (80) 

or 

Rij9 ij - - 2 H..,.l'!r = 8irGT ij9 v . (81) 
For i = j = 0,1 we can get two equations for two unknown variables N(x), n(x) 

W° - IgooRg 00 = 87rGT og 00 (82) 

»w ot ( 8 sin 2 nf\ „ T „„ ,„„. 

- N = A{ x+ ^r) N ^ (83) 

Rug 11 - ^ii-Rfl 11 = SirGTng 11 (84) 



9 „ sin nf\ ^ „„ sin nf sin nf 
x 2 + 8 5-i- C/' 2 + 2 5-i- + 4— =-f 



(85) 



where a = AirGF 2 is the coupling constant. The variation of the static energy of chiral field with 
respect to the profile f(x) leads to the field equation 

SfE s = 0, (86) 

NC (x 2 + 8^^) ,f + (x 2 + S s -^) N'Cf - (l + 42|^ + 4C/' 2 ) 7Vsin2n/+ 

+ 2(a; + 4sin2n//')iVC/' 2 + 2 ^1 + 8^^-^) (/j, — fi'x) Nf = 0. (87) 
Boundary conditions of this equation are [1, 3] 

/ (0) - 7T, (88) 

/ (oo) - 0. (89) 

To describe nucleon states, we need to quantizate classical skyrmion. The useful method that is 
used is a collective quantization. Consider collective coordinates are A (t) = ao (t) + ira(t) with 
conditions a^a^ = 1 and a^a^ = (/i = 0, 1, 2, 3 ). By this way, skyrmions to be real particles (as 
nucleon). With A(t), the U n (x,t) will trasform as 

U n (x,t)=A(t)U$(x)A- 1 (t). (90) 
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Putting Eq. (90) into Lagrangian of Einstein- Skyrme system (61), we can get 



C = -M + XTr (^4 _1 ) = -M + 2A ^ a 2 , 

[i=0 

where M — E$ is classical skyrmion's mass (in the unit c = h= 1) and 

2tt 



A = 



A = 



OO 

I NC 



Following [3], the Hamiltonian is defined as 



li=0 



Eigenvalues of the Hamiltonian (94) are 



E = M + 



1(1 + 2) 



8A ' 

where I = 21 = 2 J and (I, J) are respectively an isospin and a spin quantum number, 
nucleon (J = 1/2) and the delta (J = 3/2) we get 



M N — M + 
M A = M + 



_3_ 

8A' 

15 
8A' 



From Eqs. (73), (74), (92), (93), (96), (97), c and F v can be found as 



1 4/ A (M A - M N ) {5M N - M A ) 



n 4V6 



The baryon current is defined as 

iivpa 1 



7r 2 F 3 

167ri^r 
5M N - Ma ' 



: Tr (u- n d v u n u- n d p u n u- n d !7 u n ) . 



247T 2 n J 

From Eq. (90), the baryon currents zeroth component (baryon number density) is found 

B o_ 11 f sin 2 nf 
2tt 2 N r 2 

Using Eqs. (68), (88), (89) the baryon number becomes 



00 u 

B = J A-nr 2 ^B°dr = ~ J sm 2 nfdf = 1. 



(91) 

(92) 
(93) 

(94) 

(95) 
For the 

(96) 
(97) 

(98) 
(99) 

(100) 
(101) 
(102) 
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The isoscalar mean square radius is denned in terms of the baryon number density 

oo oo oo 

(r 2 ) /=Q = f r 2 J~gB {) ^r 2 dr = f 47rr i NB°dr = l —- [ x 2 f sin 2 nfdx. 



The isoscalar magnetic square radius is defined as 

oo 

/ r 4 ^B°4irr 2 dr 

V /M,7=0 5 



(r 2 ) I=0 



The other components of baryons current arc found 



From [3], we can get 



2tt 2 



B l = 



-g r' 



A- 1 At, 



e ijk 1 sin 2 nf t ,__ J 3 



2tt 2 \f—g r 2 
The time component of vector current V is defined as 



fxk zt- 



q i sin 2 nf 
3 n 2 r 2 



F 2 + 4 Cf' 2 + 



sin 2 nf 



Tr 



(AA-'n) , 



or 



V 



_ sin 2 nf 



3n 2 r 2 



F 2 + -1 C/' 2 + 



2 sin 2 n/ 



2 2 



2r' 



The density of charge of nucleon is defined by Gell-Mann and Nishijima 



oo 



Q = I Pnudr = h + -B, 



where 



f „ — „ /■ 2f sin 2 nf 
B = / 4:nr 2 y/^B°dr = - I — - 



/a = 



f 2 + - cr + 



sin 2 n/ 



2f dr - 



So 



47rsin nf 
3^ 



F 2 + - C/' 2 + 



2 sin 2 n/ 



n*r 



T 3 f sin 2 n/ 
2r Tr 



(103) 



(104) 

(105) 
(106) 

(107) 
(108) 

(109) 

(110) 
(111) 

(112) 



For T 3 = 1/2 and T 3 = — 1/2 we get respectively the density of charge of proton and neutron. The 
isoscalar magnetic moment and isovector magnetic moment are defined as 



-.9 r 



(V x F) 



47rr dr, 



F d 47rr 2 rfr. 



(113) 
(114) 
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According to [3], final results of Eqs. (113), (114) are defined 



1 



4 

5 Conclusions 



iM N( M.-M N) (S) I __„- 1 ^ 



(115) 
(116) 



We have constructed the term of mass of pion and calculated numerically in a case of n = 4. We 
have applied this model to define static properties of nucleon. A few obtained results are better 
than ones that predicted before [1, 3, 4, 9] as the axial coupling gA =1.22 (in [1]: gA — 1.3, in 
[3]: g A = 0.61, in [4]: g A = 0.65, in [9]: qa = 1-27), and /J, neu (mag) =-1.57 (in [1]: fjL neu {mag) = 
— 1.17, in [3]: n neu {mag) = —1.31, in [4]: /j, neu (mag) = —1.24). We have constructed successfully 
the formalism of supersymmetric composite Skyrme model. The proof was that the numerical 
solutions of field equation (24) were obtained. This was a role to calculate some static propertities 
of nucleon. The formalism of gravitational composite Skyrme model was constructed. In fact, we 
have not calculated numerically the field equation and some properties of nucleon, but we believed 
that they will be perfomed in the near future. We also believed that our results could be applied 
in some fields of cosmology. 
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figures for a case of pion's mass 




Fi gure 1: F (f) as a function of the radial distance r 
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Figure 2: Neutron and proton charge densities p(r) as functions of the distance r(fm) 
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Figure 3: The numerical solution f(x) at a = 1, /3 = 0, n = 2, 3 



